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We discuss massive scalar field with conformal coupling in Friedmann- 

Robertson- Walker (FRW) Universe of special type with constant electromag- 
netic field. Treating an external gravitational-electromagnetic background 
exactly, at first time the proper-time representations for out-in, in-in, and 
out-out scalar Green functions are explicitly constructed as proper-time in- 
tegrals over the corresponding (complex) contours. The vacuum-to- vacuum 
transition amplitudes and number of created particles are found and vacuum 
instability is discussed. The mean values of the current and energy-momentum 
tensor are evaluated, and different approximations for them are investigated. 
The back reaction of the particles created to the electromagnetic field is es- 
timated in different regimes. The connection between proper-time method 
and effective action is outlined. The effective action in scalar QED in weakly- 
curved FRW Universe (De Sitter space) with weak constant electromagnetic 
field is found as derivative expansion over curvature and electromagnetic field 
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I. INTRODUCTION 

It is quite well known fact that quantum field theory in an external background is, 
generally speaking, theory with unstable vacuum. The vacuum instability leads to many 
interesting features, among which particles creation from vacuum is one of the most beautiful 
non-perturbative phenomena in quantum field theory. Furthermore, in interacting theories 
the vacuum instability may lead to quantum processes which are prohibited if the vacuum 
is stable. One ought to say that all the above mentioned peculiarities can not be reveal 
in frames of the perturbation theory with regards to the external background, one has to 
treat it exactly. The latter has been realized long ago by Schwinger [0 on the example of 
quantum electrodynamics in the constant electric field. The particles creation in this case 
has been calculated explicitly. 

In quantum field theory with unstable vacuum it is necessary to construct different kinds 
of Green functions (GF), e.g. besides the causal GF (out-in GF) one has to use so called 
in-in GF, out-out GF, and so on []2|-|]] (for a review and technical details see H). General 
methods of such GF construction in electromagnetic (EM) background have been developed 
in . The possible generalization of the formalism to an external gravitational background 
has been given in ref. ||. Since the paper |l| it is known that causal (out-in) GF may be 
presented as a proper-time integral over a real infinite contour. At the same time, in the 
instable vacuum the in-in and out-out GF differ from the causal one. One can show |7j 
that these functions may be presented by the same proper-time integral (with the same 
integrand) but over another contours in the complex proper-time plane. The complete set of 
GF mentioned is necessary for the construction of Furry picture in interacting theories, and 
even in non-interacting cases one has to use them to define, for example, the back reaction 
of particles created and to construct different kinds of effective actions. 
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It may be likely that early Universe (EU) is filled with some type of electromagnetic 
fields. For example, recently (see |8],|9| and references therein) the possibility of existence 
and role of primordial magnetic fields in EU have been discussed. From another point 
the possibility of existence of electromagnetic field in the EU has been discussed long ago in 
[IDyir] . It has been shown there that the presence of the electrical field in the EU significally 
increases the gravitational particle creation from the vacuum. In principle, this process may 
be considered as a source for the dominant part of the Universe mass. 

Having in mind the above cosmological motivations it is getting interesting to study the 
quantum field theory in curved background with electromagnetic field (of special form to 
be able to solve the problem analytically). In the present paper we are going to consider a 
massive charged scalar field with conformal coupling in the expanding FRW Universe with 
the scale factor 0(77) (in terms of the conformal time) fi 2 (^) = b 2 rf + a 2 . Such a scale factor 
corresponds to the expanding radiation dominated FRW Universe. In terms of physical time 
t the corresponding metric may be written as following: 

ds 2 = dt 2 -Q 2 {t)(dx 2 + dy 2 + dz 2 ) , (1) 

where for small times \t\ <C a 2 /b, Q 2 (t) ~ a 2 [l + (bt/a 2 ) 2 ] and for large times \t\ 3> 
a 2 /b, £l 2 (t) ~ 2b\t\ (see [fLOR). Moreover, such FRW Universe will be filled by the con- 
stant electromagnetic field (the precise form of this field is given in the next section). Thus, 
we start from the charged scalar theory in above background. Making conformal transfor- 
mation of scalar theory (which works generally speaking for an arbitrary dimensionality of 
spacetime if correspondent conformal coupling is chosen), we remain with the theory in flat 
background but with time-dependent mass. For generality, we consider such a theory in 



d-dimensions (what may by useful also for Kaluza-Klein type theories [|12|] ) . The embedding 
of such a theory to d- dimensional constant EM field will be considered. (Of course, our 
main physical interest will be related with d — 4 where the Maxwell theory is conformally 
invariant, and one can start from the FRW Universe (|I|) with constant EM field from the 
very beginning, before the conformal transformation). In addition, in d — 4 the conformal 



coupling (i.e. the choice of scalar gravitational coupling £ = 1/6) is UV fixed point of renor- 
malization group (RG) for a variety of interacting theories (for example, for charged 

self- interacting scalar theory it is IR stable fixed point of RG |]IE||). Hence, there appears the 
additional motivation to start from conformal coupling £ = 1/6, as anyway at high energies 
(strong curvatures) [0)0] £ tends to this value. Hence, in Sect. II we solve <i-dimensional 
Klein-Gordon equation for the theory in the constant EM field and with time- dependent 
mass to get complete sets of solutions classified as particles and antiparticles at t — > ±oo. 
Using them we construct all necessary GF for the scalar field as Schwinger type proper- 
time integrals. All GF have the same integrand and differ by the contours of integration 
in the complex proper-time plane. As far as we know that it is first explicit example for 
proper-time representation for complete set of scalar GF in gravitational-electromagnetic 
background (for pure electromagnetic background it was calculated in [[?] and in pure grav- 
itational background, see for example ||14|| ). Using the exact solutions obtained we discuss 
particles creation under the combined effect of an external gravitational and electromagnetic 
fields. By means of GF we discuss the back reaction of the particles created and the effective 
action. Besides, in the Sect. IV we apply renormalization group (RG) method for calcula- 
tion of the effective action in constant curvature spacetime filled by the constant EM field. 
The typical example of such spacetime is de Sitter Universe (exponentially expanding one). 
As specific theory, we consider scalar quantum electrodynamics (QED) on above- described 
background. The RG improved effective action is calculated as the derivative expansion over 
the external fields (hence, in this Sect, we consider weakly-curved spacetime filled by weak 
EM fields). Such an effective action represents the extension of the well-known Schwinger 
effective Lagrangian for the case of curved spacetime. 

In the Conclusion the summary of the results obtained is presented. Some further 
prospects and possible applications are also discussed. 
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II. EXACT SOLUTIONS AND GREEN FUNCTIONS IN THE EXTERNAL 
TIME-DEPENDENT BACKGROUND 



In this section we will present exact solutions and GF of the charged scalar field in the 
external constant uniform electromagnetic background. In addition, the former field will be 
considered in the time-dependent mass-like potential, which effectively reproduces effects of 
a gravitational background. 

The scalar field obeys the Klein-Gordon equation: 

(P 2 -M 2 n 2 (r)))(p = (2) 

where P M = id^ — qA^ q = — |e| (for electron); x° = i], <9 = d v ; <p(x) is the scalar field; the 
dimensionality of the space is chosen to be equal to d = D + 1, d > 2 and the Minkowswski 
tensor has a form r/^ = diag (1, —1, —1, . . .) . The time-dependent potential term is chosen 
as Q 2 {?l) = b 2 r] 2 + a 2 . The external electromagnetic field will be chosen as the following: 
constant uniform electric field 

F 0D = E (Ffj, v = d^A u - d v Ap) , (3) 

and magnetic field ( for d > 3), which is defined by [d/2] — 1 invariants of the Lorentz 
transformations. This magnetic field is given by the corresponding components Hj (see for 
details @): 

[d/2]-l 

FjP = E H 3 (5 2 ^ - 5 2 J5 2 ^) . (4) 
For such a field we select the following potentials: 

A = 0, Ad = Ex , A* = A± = -Hj x i+1 6 il2j -i, J = 1, • • • , [d/2] - 1, i = 1, . . . , D - 1. (5) 

In this case solutions <p of the equation (0) are related with ones <p' with some other choice 
of the potentials A'^ for the same electromagnetic field, by the relation 

ip' = exp (iq £ (A„ - A'^) dx^j <p , (6) 
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where the integral is taken along the line: d v f* c (A^ — A'^j dx^ = A v — A' u . 

Let us discuss first the physical motivations behind the Eq. (Q). For b 2 = 0, a 2 = 1 we 
have c?-dimensional scalar field in the constant electromagnetic field. If in this case d > 4 we 
have typical Kaluza-Klein type theory in external EM background (for review of quantum 
Kaluza-Klein theories, see, for example ||12| ). 

The case d = 4 , and b 2 ^ is of special interest for us. In this case one can consider 
charged scalar field in the conformally-flat Universe (with scale factor fi(^) = y/b 2 r]' 2 + a' 2 ) 
filled by the constant EM field. Choosing the scalar-gravitational coupling constant 
£ = 1/6 (conformal coupling) and making the standard conformal transformation of the 
gravitational metric and scalar field, we come to the theory in flat spacetime with time- 
dependent mass. The corresponding field equation is given by (|2]) (EM field should not be 
transformed under the conformal transformation). Note that £ = 1/6 is UV fixed point 



of RG for a variety of theories |]T3|,|T2j, so such conformal transformation may be used also 
for interacting theories (let us remind that (y? + \p) 2 -term is conformally invariant in the La- 
grangian). Hence, d — 4 case in Eq. (Q) actually corresponds to the quantum scalar field 
in the expanding FRW-Universe with constant EM field. However, for generality we leave 
d to be an arbitrary integer number. 

As was already said, when quantum fields are considered in time-dependent backgrounds 
(electromagnetic or gravitational ones) one has to construct different Green functions. To 
this end one has to find special sets of exact solutions of the equation (||). Here we are going 
to describe such solutions. The functions <p(x) can be presented in form: 

[x\\)(p n (x ± ) , (7) 

where nonzero x l ± = x % , i = 1, . . . , D — 1, x \\— x ^-> ^ = 0,D and 

£±<Pn(?_ L )=P 2 ± <p n (x_ L ), Pl = P\ i = 1,...,D- 1, (8) 



<^* n (x_ L )f n/ (x ± )dyi ± = 5 nn , 
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(9) 



Here n = (ni, . . . , n[d/2]-i,Pi,P3, ■ ■ • ,P2[(d-i)/2]-i) is complete set of the quantum numbers 
in the space x±_ . The dimensionality of n is equal to d — 2. The expression P \ is defined 



as 



[d/2]-l 



0, d is even 



<jj ; 



Vj + Mo, ^0 
i =1 I v 2 i _ 2 , d is odd 

15^-1(2^ + 1), nj = 0,1,-.. , Hj^O 



(10) 



pl-i+pl 



■ZJ-L ' r-2] > 

where in the presence of the magnetic field some momenta p 2 j have to be replaced by the 
discrete quantum numbers rij. The number of these momenta p 2 j corresponds to the number 
of nonzero parameters Hj. Note that n includes the discrete nj and continuous P2j-i '■ 



P2j-lfn{x±) = P2j-lfn{x±) 



(11) 



Let us write 



where 



-ip D x- 



<2k 



¥p D n(x°) — Vp D n( x ° ,Pz)\p z =0 i 



(12) 



and (f(x , p z ) is a solution of equation 
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*J=) ~ (Pz - pfj) 2 ~ 



¥p D n \X°,Pz) = 



with 



(13) 



fj = x° - - 2 qEp D , p 2 = (qE) 2 + (bM) 2 , pX = (aM) 2 + p 2 D M£ + Pi . 
p z p z 

Orthonormalized and classified as particles (+) and antiparticles (-) at x° — > ±oo solutions 
of the equation ([13]) have the form |19| 



+<p PDn (x°,p z )=BD v [±(l-i)Z], ±<p PDn (x°,p z ) =BD_,_ 1 [±(l + i)e] , 

£ = ^ (P*7 - P.) ,» = | ~ \ , B = (2p)-W exp , (14) 
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18 1 for additional arguments advocating such a classification. 



One can find decomposition coefficients G (x\j °f the out-solutions in the in-solutions, 
VOO = + ¥>(z)G - -<p{x)G , C = ± • (15) 
The matrices G (cl^ ) ODe y the following relations, 

G(cl + )G(cl + ) t -G( c |-)G(c|-) t = a, 

C7 (+1+) C7 (_|+) + - ( + |-) G? (_|-) f = . (16) 

The latter can be derived from the orthonormality conditions. One can easily see that 
G M^') are diagonal, 

G (d C ')^ = V^(cl C ') » / = (PD^),/' = (P^^), (17) 

where 

#(c| C ') = C^nO^P*) « 5 f V PD n(x°,P,) • (18) 

As was already said, in case of theories with unstable vacuum one has to study different 
types of GF, which we define via A, A c , A T , A| n , Af n , and so on, according [0-0. For 
example, the field theoretical definitions of causal out-in, in-in and out-out GF have the 
form 

A c (x,x') = c~H < 0,out\T(f)(x)$(x')\0,in >, c v =< 0, out|0, in > , (19) 
A c in (x, x') = i < 0, in\T<f)(x)^{x') |0, in > , (20) 
A c out {x, x') = i < 0, out\T<f)(x)^{x') |0, out > , (21) 

where 4>(x) is a quantized scalar field, |0, in > and \0,out > are initial and final vacua and 
c v is vacuum to vacuum transition amplitude. In principle, one can calculate all of GF using 
sets of solutions (^) and some relations between the functions, 
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A c (x, x) = 9 (xo — x'o) A (x, x 1 ) — 9 (x'q — xo) A + (x, x) 
A(x,x')=i <fi(x) , <ft (x')\ = A~(x, x') + A + (x, x') , 



[X 
r+oo 



dpD + fpDn{x)g(+\ + T + l f*p D n(x') , 

-CO n 



A + (x,a;') = -i / dp* X> PBB (aO sU")- 1 -<„(ar') 



(22) 
(23) 



(24) 



A c in (x, x') = 6»(a;o - a;' )A in (x, a/) - 0(a/ o - x )A + in (x, x') , 

/+oo ^ 
dpD}2 ±<P PD n(x) ±P* PD n( X ') , (25) 
-oo n 

A c out (x ) x') = 9(x - x' )A~ ut (x, x') - 9(x' - x )A + out (x, x) , 



/+oo 
dpDY, ±l fPDn(x) ± ip* pDn (x') . 
-oo „ 



(26) 



Below we are going to do such calculations in the case under consideration. First, one can 
remark that the functions A T and Af n can be presented as follows 



±A T (x, x') = A c (x, x') ± 9(t(x - x' ))A(x, x') 



±A? n (x, x') = Al (x, x') ± 9(t(x - x' ))A(x, x') , 



±A* lt (x, x') = A c out (x, x') ± 9( T (x - x' ))A(x, x') 



At n (x,x') = A c (x,x')-A a (x,x') 



A c out (x,x') = A c (x,x')-A?(x,x') , 



A a (x,x' 



+ OO 



oo 
+oo 



dpDj2-v PD n(x) \g( 



+ Pp D n\ X 



/+oo 
-co „ 



(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
(33) 



The coefficients ([T^) do not depend on p z , thus, one can present the functions A T and A a,p 
in the convenient form 



A T ' a ' p (x,x') 



+oo r+oo 

dz 

-oo J —oo 



dpo 
2tt 



p -ipDV D ATAP „. 



Xfj, x^ , 



(34) 



Q 



A Q a ' P (v,x±,r]',x' ± ,z - z',p D 



/-tOO u _ 1 /~ \ 

^ r+oo ^ 

A Q = - i z2 dp z -ip PDn (r],x ±J z,p 



01- 



<Pp D n[V',x'x,z',Pzi ■ 



+00 



A Q = ~ i J2 dp z „ip PDn (rj,x±,z,p z 



9i-\ ) 9U\ 



^ +< ^ PD n(?7,a;±,^P^)^(+| + )"' g( 

where 



+ 1 ) VlnnVl'iX'^z'iVz, ■ 



(35) 
(36) 



(v, %±1 Z ) Pz) — VpD n 
1 



¥p D n (V, Z iPz 



^7T 



The functions Ag' a ' p have the form of the corresponding GF in EM background, where fj is 
the time, z is the coordinate along the electric field, the mass rriq = a 2 M 2 + p\, , and 
the potential of the electromagnetic field is A z = fj£. The functions 

(37) 



ip(rj, z) = / ip PDn (r], z, p z ) c npz dp z 
J —00 

with some coefficient c npz , obey the equation 

(id^ 2 - (id z - pfjf- P 2 ± -m 2 Q y )<p Q (rj,z) = . (38) 

Its solutions can be found |2(J in the form of decompositions in full sets of solutions of some 
first order equations in the light-cone variables, z^ = fj =F z, 



l(p(fj,z,p- 



1 



-1/4 f«P ( l 2 2\ -V- 

p 1 exp { — —z — z ) —1 — 



z + + z/ln(=F7r_)| 



7T_ = y/p~7I- = P- — pZ~ , ln(=P7T_) = hl(7T_) + i7T#(±7r_) , 



-<p(fj, z,p-) = 9(u_) + ip(r], z,pJ)g{ + \ + ) 



.<p(rj,z,pJ) = 0(-7T_) y>(»7,z,p_)(/(_| )* , 



(39) 



10 



and ± classification has the same meaning as one in (0). Such decompositions have the 
form of integrals over the parameters p_ . Then one can present the integrals (|35|) by means 
of those p- integrals, 



+00 



TA^=/ e(T*-)lfdp-, (40) 

n 

Taking into account that |g(-| + )| 2 = e _7rA , one gets 



Aq = -/ e(-n_) + fdp_, 

J —00 
r+00 _ 

A P o = - / 0(5r_) "/dp- • (41) 



It was shown in that the p_ integrals can always be transformed into integrals over the 
proper-time s. After the summation over n, which can be done similarly to the case d = 4 
@, one gets a result which can be written as 

±Ag = ± 9&Vo) Aq , (42) 

Aq = / .feds , f Q = f Q (fj,x ± ,ri',x' ± ,z- z',p D ,s) , (43) 

Aq = e(y ) / /q^s , e(j/ ) = sgn(y ) , (44) 
Jr c -r 2 -ri 

A% = / f Q ds + 9(z' - z) f f Q ds , (45) 

A?q= f fQds + 9(z - z') f f Q ds , (46) 

where 0(0) = 1/2. 
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1. Contours of integration ri, I^, T^, T c , T a 


/q = 


exp{- 


-iq f 

Jx' 


A id^}fg (p, T]',z- z',pd,s) fx(x±,x'±,s) 



(d-2)/2 



f±(v±,s) = c d n 



j ^sin(g£T,-s) 



m^)=^~ i/2 n 



(d-3)/2 



j ysin(gf/js) 



d is even, 



, d is odd, 



c d = (An)-* 1 ' 2 exp 
/q (v, V',z,Pd,s 



n 



(d - 4) + 2/±gF w coth(gF (// V 



2/± 



exp^-i-(»7 + »/)z-imgS 
smh(ps) I 2 v 



z 2 — (77 — i]') 2 coth(ps) 



whereas the following relations take place 

d 



ifofQ = ( m Q - P o) fQ ' 
Pq = K) 2 - (id z - pf,) 2 - Pi 



(47) 



(48) 



(49) 



(50) 



il™) fQ = iS (v- V') S{z- z')5(y L ) . 



(51) 



In accordance with fl3"4| ) one can calculate the Gaussian integrals over po and z for all 
the points s on the contours Fig.l. As a result one gets 
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where 



A c (x,x') = / f(x,x',s)ds , 
A(x,x') = e(y ) J f(x,x',s)ds , 
A a (x, x') = -A (1) (x, x') - A® ( x , x ') , 
A p (x,x') = -AW(x,x') + AW(x,x') , 

AW(a;,x / ) = ~ / f{x^, 8 )d8 , 

A( 2 )(x,a;') = / f r (x,x',s)ds , 

•/r 3 +r 2 -r a 

/ r (x,x',s) = ^=7 /(x,x',s) , 



(52) 

(53) 

(54) 

(55) 

(56) 
(57) 



a = e 



-27T/2 



(x + x' )s(bM) 2 + qEy D 



(58) 



(59) 



4s(6M) 2 cu 

and 7 a) is the incomplete gamma-function. Here 

f(x, x', s) = e iqA f\\ (x , x' , y D , s) f±(y±, s) , 

fw (xo, do, y D , s) = sinh ^ u i/2 ex P {*~y ( x ° + y D ~~ l \ ( x ° ~~ x '°^ Coth (P s ) 

-i(aM) 2 s + i^-y 2 D coth(ps) - \(bM) 2 s (x + x' f + 2qEy D (x + x' )} } , 

(bM) 2 , , x (g£) 2 , x 

= ^ coth(ps) + . (60) 

P P 

Here only A depends on the choice of the gauge for the constant field, via the integral 



A 



A' 



(61) 



which is taken along the line. To get the function / in an arbitrary gauge A', one has only 
to replace A by A' in the A. 
One can see that 

d 



i- f(x, x', s) = (M 2 fi V) " P 2 ) f(x, x', s) , 



lim^ f(x, x , s) = i5(x — x ) 



(62) 
(63) 
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Thus, f(x,x',s) is Fock-Schwinger function |l,|2l|. The contour r c — T 2 — Ti in (|53|) was 
transformed into T after the integration over and z. Then the results are consistent with 
the general expression for the commutation function obtained in |[22||. For appropriate choice 



of gauge and in d = 4 the function f(x, x', s) coincides with one from ||11|| . 

If b ^ 0, then the function f(x,x',s) has three singular points on the complex region 
—7i < Q(ps) < which are distributed at the imaginary axis: pso = 0, ps\ = —iix and 
ps2 = —ic2- The latter point is connected with zero value of the function to. We get an 
equation for ci from the condition uj = 0, 

c 2 tan(c 2 -7r/2)- (j^j =0, (64) 





Im 


s 













r 
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r i 




p 



Re s 



FIG. 2. Contours of integration T, 

where 7r/2 < c 2 < tt. The position of this point depends on the ratio qE/(bM), e.g. at 
bM/(qE) —>■ one has c 2 — ► vr and at qE/(bM) — > one has c 2 — > 7r/2. Notice, in the case 
E = it is convenient to put c 2 = 7r/2 + because of the contour T 2 must be passed above 
the singular point s 2 in the case as well. 

If b = 0, then uj = 1 and the function f(x, x', s) has only two singular points So and Si 
on the above mentioned complex region. In this case the gauge invariant function fa does 
not depend on on x + x' . In this degenerate case it follows from (0), (|^) and ([46]), 
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A a (x,x')= [ f(x,x',s)ds + 6(-y D ) [ f(x,x',s)ds, (65) 

Jv a Jr 3 +r 2 -r a 

A p (x,x')= [ f(x,x , J s)ds + 6(y D ) [ f(x,x',s)ds . (66) 

Let us return to more interesting case b ^ 0. Our aim is to demonstrate that func- 
tion x') from (|57j ) can also be presented via a proper-time integral with the kernel 
f(x, x', s) as it was done for all other A-functions. To this end let us transform the contour 
r 3 + T 2 - T a in (0) into two ones: (see FIG.2) and + T r (see FIG.3). 
Im s 



e 
e 



r 



Re s 



2p 



r 



a. 



p 



FIG. 3. Contours of integration Tjj, T; and r r 



The radius of the contour r° tends to zero. The contour T[ + T r is a infinitesimal 
radius clockwise circle around the singular point s 2 . However, it is convenient to present 
it as a combination of two semicircles Ti and T r placed on the left and the right sides of 
the imaginary axis respectively. The argument args' of the Ti radius is in the interval 
tt/2 < args' < 37r/2 and of the T r radius is in the interval — 7r/2 < args' < ir/2. Then 
can be rewritten in the form 



A (2) (x,x') 
r(x, x') = 



r,+r r 



f r (x, x', s)ds + r(x, x') 



fr(x, x', S)ds 



(67) 
(68) 
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Taking into account ( |B1| ) one gets 



7T 



f r ( x : x , s i 
\ P 



) 7Z*o const ' exp { - i^ ( x ° ~ x '^ 2 } ■ 



Hence, one can see that r(x, x') = 0, dor(x, x') = at any x — x' Q . Moreover, using (|62|). 
it is easy to see that the distribution r(x,x') obeys the equation (||). Thus, r(x,x') is equal 
to zero identically. The function f r (x,x's) is 2ir periodic function of the argument args' of 
the Ti and r r radiuses. One needs to take into account the asymptotic decomposition ( |A6| ) 
which is valid in the region — 3n/2 < arga < 3ir/2. Then, using ( |A6| ) one gets from ( |6"TD 

Jr ; +r r f(x,x's)ds , -5tt/4 < < -3tt/4, 
Jt,-iy f(x,x's)ds , -3vr/4 < (3 < -tt/4, 



A (2) (x,x') = i < 



(69) 

Jr i+ r r , -tt/4 < /3 < tt/4, 



/r,-r P f(x, x's)ds , vr/4 < (3 < 3tt/4, 

where (3 = arg (x + ig)c 2 (&M) 2 (-i) + pqEy D . 

One can verify that expression ( |69"1) is continuous in the boundaries of the f3 intervals. 
Then, using (^), one can demonstrate that the representation (^) obeys the equation 
(0). One can also verify that the representation A^'(x, x') (|56D obeys the equation (||). 



Thus, all the A-Green functions considered here, excluding those marked by the index "c" , 
are solutions of the equation (|2|). The important difference between basic Green functions 
A c (x,x r ), A« (x,x') and A(x, x'), AW (x, x') is that the first ones are symmetric under 
simultaneous change of sign in x , x' , xd, x' d and the seconds ones change sign in this 
case. 

Note finally that using proper-time kernel f(x,x',s) ( |59"D one can easily construct 
Schwinger out-in effective action 



out— in 



/poo 
dx J s _1 /(x, x, s) ds; . 



Similar out-in effective action (but in another approximation) will be discussed in Section 
IV for scalar electrodynamics. As regards to in-in effective action its representation is more 
complicated @J^] and will be discussed in the next publication. 
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III. VACUUM INSTABILITY AND BACK REACTION OF PARTICLES 

CREATED 



All the information about the processes of particles creation, annihilation, and scatter- 
ing in an external field (without radiative corrections) can be extracted from the matrices 
G (cl^') (HI)- These matrices define a canonical transformation between in and out creation 
and annihilation operators in the generalized Furry representation PH, 



a\out) = a ] (in)G ( + | + ) + b(in)G 
-b(out) = a\m)G (+|~) + b(in)G (_|~) . (70) 

Here a] (in), b](in), ai(in), bi(in) are creation and annihilation operators of in-particles and 
antiparticles respectively and a\ {out) ,bj (out) , ai(out), bi(out) are ones of out-particles and 
antiparticles, / are possible quantum numbers (in our case I = PD,n) ■ For example, the 
mean numbers of particles created (which are also equal to the numbers of pairs created) by 
the external field from the in- vacuum \0,in > with a given quantum number / is 

N t =< 0, in\a}(out)ai(out)\0,in >= g(_| + ) 2 . (71) 



(for a review of gravitational particles creation, see |24jj29[].) The standard space coordinate 
volume regularization was used to get the latter formula, so that 5(pj — v'j) ~ > & Vi , v '.- The 
probability for a vacuum to remain a vacuum is 



P v = \c v \ 2 = exp + iVi)j 



(72) 



where |0, out > is the out- vacuum. 

Remember that we are discussing the case in which the electric field acts for an infinite 
time. However, one can analyse the problem in finite times T = x° out — x° n , acting similar to 
|18| . Then the mean numbers of (pD, n )- particles created by the external field are 



N PD n = \g(-\ + )\ 2 = e-*\ if v/pT > 1, and JpT > A, and p 2 T > \qEp D \ , (73) 

where A is defined in (|13"D. The latter conditions take place for large T. 
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At d=4 ( [73|) coincides with the one obtained in ||10|| , and at b = (the gravitational field 
is absent) it coincides with the one obtained in flIB| . (bM) 2 /p 2 < 1, and the number N pon 
depends from p 2 D more weaker than from other quantum numbers ~P\. 

If the condition p 2 D (bM) 2 /p 3 << 1 takes place (the gravitational field is in a sense 
weaker than the electric one) the pr, dependence of the mean numbers ( |75| ) is similar to 
the pure electrodynamical case. Thus [pl|] one can evaluate that / dps = (eE)~ 1 p 2 T '. Then 
one can estimate the particle creation per unit of time similarly to [plj. In strong enough 
gravitational fields time dependence of the effect is nonlinear one and demands a special 
study 

To get the total number N of particles created one has to sum over the quantum numbers 
n,pE>. The sum over the momenta can be easily transformed into an integral. Thus, if b = 



one gets result presented in [18]. If b ^ and d = 4 the total number of pairs created per 



space coordinate volume has the form 



n = — , , — = — — — t exp < — 7r ) , (74) 

/ dx 8tt 2 bM F 1 p J ' v ; 



where 

qH 



(3(n) 



smh(nirqH/ p) 

The observable number density of the created pairs in the asymptotic region xq = Xg"* — > oo 
is given by 

n cr = h/Q 3 {x ) . (75) 

These results coincide with ones in [ [Hfl . 

In case 6^0 the expression ( [74] ) is growing unlimited. In this case the particles are 
created in main by the electric field, whereas the parameter b plays a role of "cut-off" 
factor, which eliminates creation of particles with extremely high momenta along the electric 
field. One can see that from the expression (|73|). Thus, the limit b — > corresponds 
to the case of the electric field which acts an infinite time when the number of particles 
created is proportional to the time of the field action. In fact, we are interesting in the 



case when the field action time T is big enough to obey the stabilization condition, which 
has the form T » (qE)^ 1 ^ 2 for intense field. As was already remarked above, in this 
case / dpo = (qE)~ 1 p 2 T. Then it is clear that parameter b has to be interpreted as a 
quantity which is inversely proportional to the field action time T making the substitution 
b- 1 -> TM(tt^p)- 1 . 

The vacuum-to- vacuum transition probability can be calculated, using formula (fT2|). 
Thus, we get an analog of the well-known Schwinger formula flj in the case under consid- 



eration. For the case 6 = the result was given in ||18|| . For the case b ^ and d = 4 one 

gets 

(aM) 2 1 



F v = exp < —fin / ax > , fi—/ J n ; lN ^ /90/lN exp 



i=0 



(Z + l) 3 / 2 /3(l) 



-/vr- 



(76) 



This result coincides with the one obtained in ITT 



Let operator of current of scalar field operator (f)(x) obeying the Klein-Gordon equation 
has a form 

1 



3n 



4>\x), <f>(x) 



(77) 



where P n = i-^ - q A^(x) and P\ 



qA^x'), and operator of Chernikov-Tagirov 



metric energy- momentum tensor (EMT) of the field operator [25| reads 



T, 



rpcan 



1 (J.V *T" "IMS J 



B 



V), 






<j>(x) 








+ 



B 



fj.ii 



"[HI 



P'\P V + P'lPp - (P'IP» - M 2 n 2 (x ) 



c, 



)1U, 



4>\x), <j>{x) 



a 



(78) 
(79) 
(80) 
(81) 
(82) 



where T™ n is canonical EMT operator. Here <f)(x) is scalar, j n is vector, and T ul/ is tensor 
under proper homogeneous Lorentz transformation. For d = 4 one can get operators of scalar 
field the current vector and EMT t^ v under general coordinate transformation by 

using formulas 
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■ flV 



(83) 



We are going to discuss the following matrix elements with these operators 



< 3n > c =< 0, out\jn\0, in> c v l 



< > c =< 0, out\T^\0, in> c v l 



< 3» > m =< 0,m|j M |0,m > , 

< > in =< Q,in\T IM/ \0,in> , 

< Jm > out =< 0, outlj^O, out > , 

< > out =< 0, out\T^\0, out > 



(84) 
(85) 
(86) 
(87) 
(88) 
(89) 



Using GF which were found before, one can present these matrix elements in the following 
form, 



<3» > c =q(P, + P'l) H)A c (x,x') 
<T^ > c = B^(-t)A c (x,x')\ x=xl + C^H)^ c (x,x) 

< 3, > in =< 3, > c + < 3, > (1) + < U > (2) , 

< 3, >° ut =< 3, > c + < 3, > (1) - < Jm > (2) , 

< > in =< T uu > c + < T uv >« + < Tuv >® , 



- /IV 



- /J.U 



< > out =< T^ v > c + < > (1) - < >W 



<3» > {1 ' 2) =q(P» + P':) H)± M 



[X, X 



(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 
(97) 



<T, V >^= B, v (-i)A^ 2 \x } x>) _ i + C flv (-i)A^ 2 \x,x) 



where GF are given by eq. ([52]), (|56|) and (0), and the relation 



A c (x,x) = - A (x, x) — A + (x, x) 



is used. 

The components < >^-' 2 > and < T^ u >( 1 ' 2 ) in the expressions ( |9"2"D and fl9~3| ) can not be 
calculated in the frame of the perturbation theory with respect to the external background 
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or in the frame of WKB method. Among them only the term < > m was calculated before 



and only in the pure electric field in flat space (6 = 0), see ||. The only expression (|91| ) for 
< > c has to be regularized and renormalized. The expression ( p0|) for term < j M > c is 
finite after the regularization lifting. The terms < j M >( 1>2 ) and < T^ v >( 1,2 ) are also finite. 
That is consistent with the fact that the ultraviolet divergences have a local nature and 
result (as in the theory without external field) from the leading local terms at s — > +0. The 
nonzero contributions to the expressions < >( 1>2 ) and < >( 1>2 ) are related to global 
features of the theory and indicate the vacuum instability. 

Let us introduce the normalized values of current and EMT (which maybe easily con- 
nected with observable values), 

j; r =]; r /n 3 (x ), (98) 
TZ = fZ/&(x ), (99) 

(in some convenient asymptotic region xq = Xq s ), where the corresponding densities per 
space-coordinates volume are 

sr = w< j» >; A - < >-> , „ = doe, 



out\ 



T^ = J ^^^ Xo = x -, (101) 



according to the definitions (|86|) - (|39"1). 

Then, using representations (^) - (|97|) one gets from (|100|) and ( |101|) , 

jZ = 2<j,>W, x = x a s , (102) 
f^ u = 2<T, u >^, x = x a s . (103) 

To study the backreaction of particles created on the electromagnetic field and metrics 
one needs the expressions < > m and < > m at all times Xq. Below we are going to 
analyse these expressions for different times. Note, that the functions A c (x,x') (|52| ) and 
A^>(x,x') Q55|) at x = x' are even functions on xq. Thus, the functions < T^ v > c and 



< T^y >^ are also even ones and do not vanish at xq — > 0, whereas the functions < j M > c 
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and < >W are odd ones and vanish in this limit. Moreover, we have < > c = for all 
x at b = 0. 

The proper-time integral A^ 2 ^(a;, x') ( 57|) is a odd function on xq at x = x' and vanishes 
at xq — » 0. Thus, the expression < T^ u >^ is also a odd function on x and vanishes in this 
limit. The term < j M >( 2 ) is an even function on Xq and is different from zero in this limit 
if E ^ 0. 

One can see, using the expressions obtained for GF, that all off-diagonal matrix elements 
of EMT are equal to zero: 



< T. v > c =< T, u >«=< > (2) = 0, fi + v- 



(104) 



and only x D components (along the electric field) of the currents are different from zero and 
vanish in the absence of the electric field. 

Below we are going to analyse contributions to the physical quantities (|102|) and (|103|) 
for d — 4 at xq = Xq s , which are stipulated by the GF A^ 2 ^(a;, x'). 

Using the asymptotic form of A^ 2 )(x,x') given in (|^8|) , one can see that at Xq >> 
yfp/ (bM) the following asymptotic expression takes place 



<j» >^=q 2 \E\p- l h-5l 



(105) 



Taking into account eq. ( |A5| ) one gets an asymptotic expression for diagonal < T^ v 
< Too > (2) = \p 2 x 2 + a 2 M 2 + qH coth (nqH/p) 



n 

px 



< T n >^=< T 22 > 



(2)_ 



qH coth (nqH/ p) 



n 



2px 



< T 33 



(qEx ) 2 + 



n 

px 



(106) 



2tt(6M) 2 _ 

where n cr is given by (|74]) . 

Doubling the expression (|105|) and (|106|) according the eq. (|102|) and ( |103|) , one gets the 
mean densities for current and EMT of particles created. It turns out that these quantities 
are proportional to the density of total number of particles and antiparticles created 2h cr 
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for the infinite time and do not change their structure with increasing of Xq. The latter 
means that one can consider all the expressions obtained at any fixed Xq if Xo S> ■ s fpjbM. 
For a strong background a 2 M 2 /p < 1 and therefore this time has to obey the condition 
x » a/b. Thus, in the strong background our asymptotic conformal time Xq corresponds 
to the large cosmological time t. 

Note, that one can neglect the second term in the brackets of the expression ( |106| ) for 
< T 33 >( 2 ) at bM/{qE) < 1. Also one can neglect both the term a 2 M 2 in the brackets of the 
expression ( |106| ) for < Too > in case of strong external background a 2 M 2 /p < 1 and third 
term in the same expression if the magnetic field strength is not big enough qH/p < 1. The 
current density j?f = 2 < >^ does not depend of the asymptotic time. At b — > this 
expression coincides with one for flat space < > cr = 2\q\n cr 5^. The pressure component 
along the electric field direction = 2 < T 33 >( 2 ) is growing with time upon the action 
of the field. However, if qE/(bM) << 1 then the asymptotic condition for xo is consistent 
with the fact that the term (qExo) 2 in the expressions < T 33 >( 2 ) from ( |106[ ) will not be 
dominant until big enough time xq. In this case one can neglect the contribution which 
depends on the electric field if the field is switched off before. Note, that more accurate 
analysis of back-reaction of created particles to current requires the numerical estimations 
(compare with purely EM case, |28| ). 

The components of the pressure in the directions perpendicular to the electric field 
T"-q = 2 < T11 >^ and T 2 C £ = 2 < T22 >^ are growing in relation to the T 33 when 
the magnetic field is increasing. However, a very strong magnetic field decreases all the 
components of T^ u and because the particles creation will be decreased. 

One can remark that the asymptotic behaviour of < j M >W and < T^ v >W is defined 
by the asymptotic expression for A^(x, x') (|A11|). Then one gets 



< U > (1) =< U > (2) » < T ^ > (1) =< T ^ > (2) • ( 107 ) 

The expression ( PD| ) does not need to be renormalized, thus, one can easily verify (using 
( AT| ) and ( |A2j) ) that the relation < > c ~ Xq 1 — > holds asymptotically. An estimation of 
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the finite part of < T^ IV > c can be made only after renormalization. We are going to consider 
this problem, together with others, related to the renormalization, in the next paper. 

To get an idea about the behaviour of the expressions (p^ ) and ( p4|) at finite time let us 
estimate their components for some small xq, namely for Xq << p/(bM) 2 . Dislocating the 
contour T 2 to the real axis and neglecting the divergent terms in < > c + < T^ u >^> 
which do not depend on the background, one can see that 

< j, > m = 3ft < j, > c + < j, + < j, >^ , (108) 
<T llI/ > m =^<T llI/ > ren + <T flu >^ + <T^>^ . (109) 

Here < T^ u > ren is obtained from < T^ v > c as a result of such a procedure and 

< J, > {3) = q (P, + P'l) (-t)A^\x,x')\ x=xi , (110) 
<T^ >^= B^(-i)^\x,x')\ x=xi + C^(-i)^\x,x) , (111) 

where A^ 3 '(x, x') is defined in (|A9|) . As was already remarked, < j M > c and < >^ are 
odd functions on x - They vanish at Xq — > 0. That is why the leading term in ( |108| ) is 
defined by < >^ . Using the expressions (|B~4T), obtained in the Appendix B, one gets 



< 



U > (2) = frP% (112) 



If bM/(qE) > 1, then the odd function < T M „ >( 2 ) vanishes at x — > and the leading 
contribution to dlQ9|) is defined by 9ft < > ren and < >< 3 > . However, HbM/(qE) « 1, 
in the domain where qE/(bM) 2 » Xq » (qE)~ x , the situation is different and the leading 
contributions can come from < T^ u >^ . That happens since the expressions 9ft < T^ u > ren 
and < T^ v >® remain finite at b — > 0. Using ( |B4| ) and ( |B5| ) one can find expression for 
diagonal < TL, > (2) , 



<T 00 >^= £ < T U > (2) WM'xorr ( 

< T n >( 2 )=< T 22 >( 2 )= qHcoth(iTH/E)x 7t n cr , 

< T 33 > (2) = 5qEx h cr . (113) 
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^From it one can see that < T 33 >( 2 ) and < T 00 >^ are growing unlimited at b — > because 
of increasing of n cr . 

The expression ( |112|) at bM/(qE) << 1 coincides with the asymptotic one (|105|) . That 



demonstrates that in the quasi-flat metrics, when the particles are created mainly due to 
the electric field only, any time far enough (x + T/2 >> (qE)^ 1 ^ 2 ) from the initial time 
(x™ = —T/2), can be considered as the asymptotic time. However, that is not true for the 
EMT < >W and the expression ( |113| ) differs essentially from the asymptotic one ( |106| ) . 



That happens since the vacuum definition for the particles with big longitudinal momenta 
G°3 > qExq), differs essentially from |0, out > vacuum. Namely, those particles mainly 
contribute to < >^ . Technically that means that one can not obtain the normal form 
by means of the A ou *-function and, consequently, to calculate the mean value of EMT of 
particles created at xq « \fqE jibM) it is not enough to know components of < T^ v > rn . 
Note finally that the expressions for mean values of energy momentum tensor may be applied 
for the estimations of back-reaction of created particles to the gravitational background like 
it has been done in ||26|| . However, such study involves the renormalization of EMT. It also 
cannot be done analytically. We will discuss these questions in an another place. 

IV. EFFECTIVE ACTION IN SCALAR QED IN FRW UNIVERSE WITH 
CONSTANT ELECTROMAGNETIC FIELD 

In this section (which is somehow outside of above discussion) we will consider d = 4 
scalar QED in curved spacetime with EM field. That is more complicated, interacting theory. 
Using such a theory as an example we will briefly discuss how one can generalize the results 
of above discussion to interacting theories. We again consider gravitational-electromagnetic 
background of weakly curved constant curvature spacetime. The Lagrangian of the theory 
may be written as following: 

L = L m + L ext} (114) 
L m = \ (d^ - qAM 2 + i (d^ 2 + qA^if - X -M 2 <f + ^Rcf - ^A0 4 - \f, v F^ , 
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L cx t = A + kR + aiR 2 + a-zC^vap + ctzG , 

where <p 2 = (f)^ = <p\ + </>|, the parameterization of fields in terms of two real scalars is 
taken, R is curvature, C^ap is Weyl tensor, G is Gauss-Bonnet combination. The necessity 
of introduction of L ext is dictated by the multiplicative renormalizability of the theory (see 



T2|j). Note also that we consider the theory in curved spacetime with EM field, hence 
An —>■ + A^, where A^ is background EM field. For simplicity, we limit below only in 
constant curvature space. The well known example of such spaces is De Sitter space which 
is often used for description of the inflationary Universe (exponentially expanding one). 

There are few ways to study such a theory on the quantum level. When one treats the 
external background exactly, like it has been done in previous section, the first step is to 
calculate scalar Green functions. 

Due to the fact that the Maxwell theory is conformally invariant, Green functions for 
EM field will be the same as in the flat space. Then, using the proper-time representation 
for the Green functions, quantum corrections to r out _i n or r in _ in can be calculated (where 
rout-in; Tin-in are effective actions for the probability amplitudes and for the mean values 
respectively). Then, of course, the external background is kept to be exact and perturbation 
theory over only coupling constants A, e is used. However, such calculation is extremely 
hard what can be already understood from very complicated form of the Green functions 
in the previous section. Indeed, these Green functions should be used instead of standard 
momentum space propagators in Feynman diagrams (for an explicit example in pure constant 
EM background, see f|). Hence, such calculation using Green functions of the previous 
section will be given in another place. 

Instead, we will adopt quasi-local expansion for the effective action here. We will make 
the explicit use of RG ]12| in such calculation. Note that in quasi-local approach to effective 
action (i.e. when we do not treat an external background exactly but make the derivative 
expansion of the effective action over derivatives from metric and EM field) the difference 
between r out _i n , and r in _ in is not seen. That is why we denote the effective action here 
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as T simply. (Of course, particle creation phenomenon is also hard to see in quasi-local 
approximation.) As background we will consider De Sitter space with background EM field 
(no background scalar field!). 

Using the fact that the theory with the Lagrangian (|115 ) is multiplicatively renormal- 
izable, one can write RG equation for the effective action on the above background as 
following: 



d d ~ d \ 



a A, 



(115) 



where A^ is background EM field chosen so that F^ u is (almost) constant (the well-known 
example is constant electric field), Aj = (q, A, M, £, A, k, ai, a 2 , a 3 ), 7^ is 7-function for A^, 
and Pi is /3-f unction for A«. 

The solution of RG equation ( |115| ) can be easily found using the method of characteris- 
tics (see |T5| for flat space and JTB| for curved spacetime). In such a way, one obtains RG 
improved effective action (which makes summation of all leading logarithms of the pertur- 



bation theory). Dropping the details which are very similar to the ones given in ref. |TB 



we 



get RG improved effective action up to the terms of second order on the curvature invariants: 



T RG = A(t) + K{t) + ai (t)R A + a 3 (t)G - - 



1 q 1 



F 2 

4g 2 (t) 



(116) 



where the effective coupling constants are given as following (see |15] for the flat space and 
T6|j for curved spacetime): 



q\t) = q 2 ^l- 

A(f) = Y/{t) 



C = arctan 



2q 2 t 
3(4vr)< 



$ 2 (t) = $ 2 1 



2qH 
3(4?r) : 



1 



V719 tan(-V719 logg 2 (t) +C) +19 



10A 



V719 V q 2 



19 



-V719 logg 2 



M 2 (t) = M 2 



q 2 (t) 1 _26/5 cos 2 / 5 (I V719 log q 2 + C 



cos 2 / 5 (|V719 logg 2 (t) +Cy 
q 2 (t)V 26/r ° cos 2 / 5 ^V719\ogq 2 + C 



cos 2 / 5 (f\/719 logg 2 (t) + C)' 
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7t 8t 

fl2(t) = ° 2 + 60(4^' ° 3(t) = % - 45(4^' 

X(t) = A + M^A^t), K{t) = K + 2M 2 (i- l/Q)A l (t), ai (t) = a t + (£ - l/6) 2 Ai(*) 
* dt 



~a 2 (t)V 52/5 cos 4 / 5 (iv 7 ^ logg 2 + C 

A 1 (t) - 1 



(4tt) 



q 2 



(117) 



cos 4 / 5 (|V719 logg 2 (t) + C) 

Note that one should use the fact that = 0, G = -y- on de Sitter background. The 

classical Lagrangian is used as boundary condition in the derivation of Trc- 

Now, the question is about the choice of RG parameter t. We are dealing with the theory 
with few effective masses. RG improvement in such a theory is not easy due to the fact that 
one has generally speaking to generalize the total mass matrix (there are two masses from 
the scalar sector plus four more from EM sector, see |16|jl7| ). 

In order to find RG parameters t one has to specify situation in more details. For 
example, let us consider the case when the background EM field is chosen to be constant 
magnetic field, i.e. j\F^ u F^ u = \H 2 . Then, one can consider different limiting cases. 

First, let H 3> R. Then all effective masses are becoming proportional and there is 
unique choice for RG parameter t : t = |ln^. In this case, the curvature effects are not 
relevant, Trg is given by its last term in ( |117| ). That is the case actually discussed by 
Schwinger jlj]. 

Second, let R^> H. Then, the unique choice for t is t — ~ In Now, EM field effects 
are not relevant. We get purely gravitational effective action discussed in ref. |12|JT3|JT6| . 



When both fields are of the same order the choice for t (in order to make summation of 
leading logarithms) is t = \ In R ^ A ^ 2 qH ■ If, in addition, we choose the initial values for mass 
and £ as following: M 2 = 0, £ = |, then we get 

St 2 / 2qH \ H 2 

= 'm(^Y 2 ~V~ J ~z ' (n8) 

where only H-dependent terms kept. The expression (ITS) generalizes the well-known ef- 
fective Lagrangian for magnetic field obtained by Schwinger [[]J to curved spacetime. Note 
that in ( |118| ) the curvature effects have been taken into account. 
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It is known that in flat space |T[ the Schwinger effective Lagrangian leads to the stationary 
point H ^ (due to quantum corrections). However, this stationary point is known to be 
the maximum of the effective action. Hence, it does not lead to the possibility of new 
ground state with account of quantum corrections. In curved space, one can analyze Trq 
( |118| ) using equation of motion . The result of this analysis shows that position of the 
flat-space maximum H ^ is moving due to curvature corrections. However, the nature 
of this stationary point is the same (it is the maximum of the effective action). Similarly, 
one can analyze (now numerically) the case with M 2 ^ and an arbitrary £. However, we 
expect that existence of new ground state may occur only in the regime with strong curvature 
(where an external gravitational field is treated exactly like in the previous section). 



V. CONCLUSION 

In summary, we considered massive scalar field in expanding radiation dominated FRW 
Universe filled by the constant EM field. Taking scalar-gravitational coupling constant 
to be equal to its conformal value and making conformal transformation, the theory is 
formulated as flat space theory with time-dependent mass and in external EM field. For 
such a theory, which is generally speaking theory with unstable vacuum, we found proper- 
time representation for all Green functions (i.e. for in- in, out-in, and so on Green functions). 
Note that proper-time representation (over a finite contour in proper-time complex plane) 
for in-in Green function in the combination of the external gravitational and EM fields is 
given here for the first time. 

As some applications, we discuss particles creation by the external fields in arbitrary 
dimensions. Combined action of gravitational and EM fields may produce new interesting 
features in this phenomenon, in particular, affect the rate of particles creation. 

Using the proper-time representation for Green functions, the proper-time representation 
for the effective actions is also found (or, in other words, vacuum polarization is found) 
keeping external background exactly. 
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Let us discuss one more possible application of the Green functions obtained-the calcu- 
lation of radiative corrections to the transition amplitudes and to mean values of different 
quantities in the interacting theory. The corresponding general theory (Furry picture for 
theories with unstable vacuum) was developed in [|3|-[(J . As an example, let us consider scalar 
self-interacting theory with the interaction — ^y(0 + 0) 2 . The same background as in Sect. II 
will be considered. Let us imagine that one has to calculate first-order radiative corrections 
to the two-point Green function which serves for calculation of mean values. Schematically, 
we have to consider the following diagram; 




FIG. 4. A radiative correction to the propagator 

with Ai n functions. Hence, these corrections are proportional to in-in propagator in co- 
inciding points, i.e. ~ Af n (x,x). Or more completely, above-drawn diagram is proportional 
to the expression 

~ J A{y,x')A{x',x )A{x' } z)dx. 

The divergent part of the diagram is well-known, it is the same as in flat space and is 
proportional (in dimensional regularization) to . The finite part of the diagram can 
be easily found, using explicit form of the proper-time representation for the corresponding 
Green function, given in Sect. II. Of course, this finite part is different from one in the 
absence of the external background. To calculate the finite radiative corrections to the 
causal propagator, one has to take the same diagram with A c functions (or with out-in 
Green functions). The divergent part of the diagram will be the same, however, the finite 
part is different from the previous case (as it follows from the explicit structure of Green 
functions, presented in Sect. II). In the same way one can analyse other radiative corrections. 

Using explicit form of Green functions in Sect. II, we calculate the effective actions 
(vacuum polarization) in proper-time representation. The knowledge of out-in effective 
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action gives an alternative way to define the particles creation a la' Schwinger [jlj (for an 



explicit example see ||12|| ). The in- in effective action can be used to study the back reaction 
of the particles created to the external background. Such an analysis is not easy and will be 
presented in another place where also a generalization of results of Sect. II for an arbitrary 
£ will be done. For example, using explicit form of in-in GF in proper-time representation 
it could be of interest also to construct proper-time representation for in-in effective action. 

In Sect. IV we, presented another approach to the effective action (derivative expansion 
of effective action) in the external gravitational- EM background. Scalar QED is considered 
as an example; RG improved effective action (up to the terms of second order on curvature 
and EM strength) is calculated on constant curvature weakly curved spacetime with weak 
constant EM field. Such an effective action gives the extension of the well-known Schwinger 
effective Lagrangian, taking into account curvature effects. It may be also applied to the 
study of back reaction of quantum field theory to external background. 

Finally, similar technique may be applied to analyze the behaviour of spinor fields in 
gravitational-EM background. The calculation of all the Green functions in such a theory, 
using proper-time representation, may be the necessary step in the study of chiral symmetry 
breaking in QED and in the four-fermion models under the action of gravitational and EM 
fields. Such a study may have an immediate important application to early Universe, for 
example, through the construction of inflationary Universe where role of inflaton is played 
by the condensate < >. One can also analyse symmetry breaking phenomenon under 
the combined action of gravitational and EM fields in the Standard Model (using also its 
gauged NJL form ||27||), or Grand Unified Theories in the same way as it has been done in 
curved spacetime (without EM field) 

Note also that GF investigation developed in this paper maybe extremly useful for the 
study of Casimir effect due to combined action of gravitational and electromagnetic fields 
(for an introduction to Casimir effect in pure EM or pure gravitational case, see |30|j). 
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APPENDIX A: ASYMPTOTICS OF A- FUNCTIONS 



Let us calculate the asymptotic behavior of A^'(x, x') and A^ 2 ^(x, x') in the case Xq >> 
p/ (bM) 2 , and x — > x' at d = 4. 

First, it is useful to take into account the following formulas 



Pi f{x, x', s) = P>*f(x, x', s), 3 = 1, 2, 3 



P P' f(x,x's) 



Pj + 2 lu- l {bMysx ) -iuj- L (bMys 



f(x,x',s) 



(Al) 
(A2) 
(A3) 



Then, using eq.(|X3|) and equation (|5"2"D one gets 



P P'lf(x,x',s) 



d 



i— + Pf + M 2 n 2 (x ) + 2 (^(bMfsxv) - uj~H(bMf 



f(x,x',s) 



(A4) 



If the A-function obeys the equation (|2|) then the action of the operator — ij- on this function 
is equal to zero. Further we are going to use the formula (|A4|) for simplification of the 
calculation. Besides, one can see that for such kind of functions the following relation holds, 



(>*P> _ M 2 Q 2 {x )) A ( -\x,x')\ x=xi = 2d 2 A(-\x,x) . 



(A5) 



So, if one uses these formulas in calculations it is enough to find GF asymptotics for the 
case xo — x' Q = 0. 

Let us find the asymptotic behavior of A^ 2 ^(x, x') function given by the eq. (|5"TD. lib ^ 
and Xq — x'q = the kernel f r (x,x's) has no singular point Si = —in /p. Below the line of 
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contours 1^ — T a next singular point of this function is S2 — in /p. Thus, one can shift the 
line of the composite contour T 3 — T a below along the imaginary axis until the neighborhood 
of the point S2 — in / p. The contour which is obtained in such a way from + 1^ — T a can 
be closed on the parts 9fts — > ±00 and then can be transformed into a closed contour which 
includes the points s\ and S2 as well. At the same time this contour can be situated far 
enough from the point s\, so that Ipsa;" 1 1 is always not zero at it. Then one can use the 
asymptotic decomposition [p3 



-f{l/2,a) = ^-e- a a- 1/2 \l + 0{a- 1 )] , x Q > . (A6) 



The contribution from the first term of eq. flAq) at the contour in question is expo- 
nentially small since 9ft (— ipsu~ l ) < 0. The rest terms of eq. ( |A6[ ) form a series in in- 
verse powers of Xq(6M) 2 / p, the corresponding functions in the decomposition coefficients 
f(x, x', s)e~ a a~ l l 2 [1 + 0(a -1 )] have only one singular point s\ (the pole) inside the con- 
tour. The contributions from these terms can be estimated tightening the contour to the 
point s\. The first one of these terms defines the leading contribution into the asymptotics 
of AW(x,x'), 

A {2 \x,x r ) = — ^= / f(x,x',s)e- a a- 1/2 ds , (A7) 

2 V 7r Jr R 

where the contour T l R (see FIG. 3) is a circle with infinitesimal radius around the singular 
point s±. Calculating the residue, one gets an expression for A^(x, x') which defines the 
leading asymptotics in < >^ and < >( 2 \ 

A(2) ^ x ') = l p { x n + x' ) exv{iqA + b E{xo + x ' o)y3 ~ Mmy + \ y±qF cot {7xqF,p) y±} 

(A8) 

where h cr is defined in (0). This expression is also valid at x < 0, since the function 
A^ 2 \x, x') is an odd one in xq at x = x' . 

Consider the asymptotic behavior of A^(x,x') function given by eq. floPP- Since pstu^ 1 
is not equal to zero and dt(— ipsuj^ 1 ) < at the contour T2 the corresponding asymptotic 
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contribution in integral fl56|) is exponentially small. Then one needs to evaluate only the 
part of this integral given by form 

A (3 W) = "^ / f(x,x',s)ds , (A9) 

Let us introduce the variable r, ps = —in + r. Since ipsu^ 1 ) < at the contours T 3 and 
r a , and since \uj\~ 1 increases monotonous with |r|, then the leading asymptotic contribution 
is defined by the behaviour of the function f(x, x', s) at small r and has the form 

A (3 \x,x') =e- l7r/i 7i- 1/2 f(x,x', Sl ) 1™ drT- x e~ lTpx l . (A10) 

Jo 

Calculating the integral one gets 

A (1) (x,x') = sign(x )A (2) (x,x') , (All) 
where A^(x,x') is given eq.( 



APPENDIX B: SMALL TIME EXPANSION OF A( 2 )-FUNCTION 

Let calculate a small time expansion of A^ 2 ^(x, x') function given by eq. ([57]) in a case 
Xq << p/(bM) 2 and x — > x' at d = 4. The small a expansion of the incomplete 7-function 



is valid at the contours r a , T 2 and T 3 and it has a form [23 



7 (1/2, a) = e - a a 1/2 [2 + (4/3)a + 0(a 2 )] , (Bl) 

A * A ~ * ~ 

where second term in the square brackets is necessary for the calculations of P'qPq and P'^Ps 
actions. Since term (|68 ) is equal to zero it is convenient to calculate the function A^ 2 '(x, x'), 
closing the contour r 3 + T2 — T a on the area — > ±00 and then tightening it to the singular 
point s 2 . Then the leading contributions to < j M >( 2 ) and < T^ v >^ 2 \ are defined by the 
integral 



A®(x, x') = — -= / f(x, x', s)e~ a a 1/2 [2 + (4/3) a] ds , (B2) 
2a/7t Jr;+r r 
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where the contour Ti + T r (see FIG. 3) is a infinitesimal radius clockwise circle around the 
singular point s 2 - According to eq. uj = at s = s 2 . Then in the neighborhood of this 
singular point if s = s 2 + s' one gets expansion 



lo = Jps' + (1/2 W') , 

2 



(B3) 



UJ 



10 



. 1 


(t 



















2 fg£y f g^V 



i + 



Um 



C2 + 



\bM ) 



Calculating residue at the point s 2 one finds the small time expansion of A (^-function, 
A®(x, x') = e iqA { [i(x + x' )c 2 {bM) 2 / p - qEy 3 ] n®/(qE)(p Q 



+ 1 



— (Xq + X c 



6qEp 



; K(2) + (l/2)(x + x' )(y 3 ) 2 qEK(0) 



(B4) 



y 



<Po = exp ( i V( x o + ^o)z/ 3 - (x - x' f 



P 3 (y 3 ) 2 

4c 2 (6M) 5 



r 2 



\a 2 M 2 /p - c-\l - 1/2) - c- 1 f 



+ (qH/p) coth(c 2 qH/p) + 2c 2 1 



L ~ \bM J 


2" 




=1 + (ffi) 


4" 




2 + (f§) 


4 



P (2) 

— — rr ; 



n 



(2) 



q 3 HE 2 p 3 ' 2 



-c 2 a 2 M 2 /p 



7T 3 / 2 C 2 



C 2 + (w) + (fH) 



(6M) 3 sm(c 2 qH/p) 



If bM/qE « 1, coefficients i^(Z) and n^ 2 ^ from (|B4j ) have more simple form, and in the case 
of an intensive electric field (a 2 M 2 /(qE) < 1, \H/E\ < 1) one has 



K(l) = -(-i)'(l/2)7r'n 



(B5) 
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